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Abstract 

We calculate the hadronic tensor for inclusive semileptonic B — > X^^v decay to 0(as). This 
allows Oiois AqcD/^ife) corrections to hadronic invariant mass observables to be directly evaluated 
with experimentally required cuts on phase space. Several moments of phenomenological interest 
are presented to order ©(a^ AqcD/w-b) and C'^AQ(jj-)/m^^ , allowing a consistent extraction of the 
HQET parameters up to O^Aq^j-j/m^^ and the h quark mass with theoretical error ~ 50MeV. 
The hadronic invariant mass spectrum is examined with a general moment to obtain observables 
that test the theoretical error estimate assigned to these parameters; in particular, fractional mo- 
ments that directly test the OPE for inconsistencies in the hadronic invariant mass spectrum 
are reported. The m;, Aqcd/?71c expansion present for fractional moments of the hadronic invari- 
ant mass spectrum is discussed and shown to introduce a numerically suppressed uncertainty of 
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I. INTRODUCTION 



Inclusive semileptonic B ^ Xciiy decay offers an opportunity to measure the Cabbibo- 
Kobavashi-Maskawa (CKM) parameter \Vcb\ and the bottom quark mass 
0, y, S]. Measurements of these parameters are crucial to the B factory program of over- 
constraining the CKM sector of the standard model Experimental studies of moments 
of the differential decay spectrum of B ^ X^i combined with a measurement of the total 
inclusive decay rate are useful in extracting these parameters, as these observables can be 
measured cleanly by experiment, and calculated from QCD without model dependence using 
an operator product expansion (OPE). 

The OPE is an expansion in powers of the ratio AqcB/^b, where the terms in this expan- 
sion are parameterized using heavy quark effective theory (HQET). The OPE demonstrates 
that in the oo limit, inclusive B meson decay spectra are equal to b quark decay 

spectra. To extract and \Vcb\ from the inclusive decay spectrum with high precision, 
one needs to accurately know the relevant matrix elements of terms in the OPE that the 
spectrum depends upon. Extensive theoretical effort has been devoted to calculating the 
decay rates and moments of various spectra; to test HQET by the extractions of these non- 
perturbative parameters from different spectra, and to obtain \Vcb\ and a precise value for 
the b quark mass. 

Experimental results have been reported by the BABAR, CLEO and DELPHI collabora- 
tions measuring various B meson decay spectra and moments H, 12, 12,^, Ih, 3- Recent 
analysis of this data 18| finds \Vcb\ = (40.8 ± 0.9) x lO^^ and m^^ = 4.74 ± 0.10 GeV, 



where the experimental uncertainties dominate the extraction. As the experimental errors 
are expected to decrease in the near future, it is appropriate to reexamine the theoretical 
error assigned in this extraction. The largest contributions to the theoretical uncertainties 
introduced in these fits come from the estimated size of the 0(as Aqcd/"^6) corrections with 
a lepton energy cut, the Aqcp/m^^ terms, and the 0(^Aqqj) /mlmfj terms introduced 
due to the HQET expansion employed for rric jl9l |. 

In past calculations, the lepton energy cut dependence of the 0{asAQCB/fnb) terms was 
not calculated, and these terms were treated as a source of error in the determination of 
nib and \Vcb\- In this paper, we improve upon past results by calculating the lepton energy 
cut dependence of the O^asAQCD/rrib) terms. With the calculation of these terms and the 
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moments presented in this paper, global fits will allow precise determinations of \Vcb\ and 
rrih to occur from the inclusive decay spectrum. 

As the precision of determinations of rrih and \ Vcb\ improves, it becomes important to test 
the consistency of the OPE more precisely. Observables that do not depend strongly on the 
nonperturbative parameters that introduce the dominant uncertainty in extractions of rri}, 
and \ Vcb\ allow one to test if the uncertainty assigned for all higher order terms in the OPE is 
sufficiently large. By examining a general moment of the lepton spectrum [2^, observables 
of this type, called OPE testing moments, have been found. In this paper, we apply this 
technique to hadronic invariant mass moments. By testing the error assigned to higher 
order effects experimentally we improve the confidence in the theoretical error assigned due 
to these effects in determinations of mi, and \Vcb\ from moments of semileptonic inclusive b 
decay. This allows extractions of from these decay to occur in a relatively theoretically 



clean and unambiguous fashion 



2l| . These results can be combined with the results for the 



lepton energy spectrum for cross checks and fits to determine the HQET parameters. 

The structure of this paper is as follows. In Section |n] the 0{as) contribution to the 
hadronic tensor is presented. Section ITTTl reports on moments in the IS mass scheme 0,13, 



and discusses the nib Aqqd /rric expansion present in fractional hadronic invariant 
mass moments. The decay width to 0{as Aqcn/i^b) and O (jiq^j^ / mfj and the error that 
should be assigned in the fit of the moments presented in this paper is discussed. The 
dominant parameters affecting the extraction of \Vcb\ inclusively, ml^ and Ai, are extracted 
from known moments. Observables appropriate to precisely test the consistency of the OPE 
are reported and moments that allow a measurement of the b quark mass with minimal 
theoretical error due to unknown matrix elements are presented. 



II. 0{as) CONTRIBUTION TO DECAY SPECTRUM 



A. Hadron Tensor Decomposition 



The 0{as) corrections to semilep 



;onic B ^ X^if decay have been known for particular 



spectra and moments for some time 2^ 



221, 



28| . The decomposition of the triple differential 
decay spectrum in terms of structure functions has not appeared in the literature to date, 
although the limit of this spectrum appropriate for a massless final state is known [2^. The 
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triple differential decay spectrum must be known to allow for the experimentally required 
cuts on the kinematic variables to be imposed in calculating the 0(q;s AQCD/^T^b) terms and 
to perform a general moment analysis, and so we present it here. 

We decompose the triple differential decay spectrum in terms of the invariant mass of the 
W boson y = /rrib^ where is the momentum of the lepton pair, the c quark jet invariant 
mass z — {rribV — qf /rrib'^ and the charged lepton energy — Ei/rrib- This spectrum is 
written in terms of a lepton tensor L^i, and the hadron tensor W^^ , 

1 dV 



To dy dz dE^ 

where 



= W^%yrz)L^.{y.z,Et), (1) 



192^^^ ■ 

Integrating over the charged lepton energy the differential decay spectrum becomes 

^ l2EotW^Ay,z)L^''{y,z) , (3) 



Fq dy dz 



where Eq = \/2{\ + z — y) is the leading order energy of the c quark jet, p = ml/ml and 



t = y 1 — z / Eq^ is the rapidity of the c quark. 

The hadron tensor can be decomposed in terms of the initial B meson four momentum Q'* 
and the hadronic decay products four momentum — — q^, where q^ is the momentum 
of the lepton pair. This tensor can be calculated from the discontinuity of the time ordered 
product of the current = C7^(l — 75)^, 



W^"" = ^Im 



: j d^xe'^^-"''''>^B\T[j'"'{x),r{0)]\B) . (4) 



This tensor is calculated by considering the quark-gluon level processes involved in this 
decay. The spectra obtained from the parton level discontinuity are expected to accurately 
describe physical B meson decay spectra so long as observables are sufficiently inclusive. 

The tensor decomposition in terms of the four vectors Q'^ and P'* yields five non trivial 
structure functions Wi, 



W^iy, z) = Wi{y, z) [P^'Q" + P'Q" - P ■ Qg^" + te^'^'^Qc^Pp) (5) 
- W^2(y, z)g^''' + W^{y, z)Q^Q'' + W^, z){P^Q'' + P'Q^) + Wr,{y, z)P^P\ 
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FIG. 1: The one loop forward scattering diagrams. The hadron tensor is derived by calculating 
the imaginary part of the diagrams. 

The operator product expansion of the structure functions is known to C(^AQ(-;j)/m^^ jl|, 
. There are two nonperturbative parameters at O^^Aq^^j^/m^^ labelled Ai^2 and six 
parameters at order O^AQQ^y/^b) ^ labelled ti^2,3,4 and pi,2 the definitions of which can be 
found in 



B. 0{as) Contributions to Hadron Tensor Structure Functions 

The hadronic structure functions Wi in Eq. (0) have the perturbative expansion: 



{y, z) = {y, z) + -^W} (y, z) + ^(a,^) . (6) 

The 0{as) contributions to the structure functions are calculated by taking all cuts across 
all intermediate state contributions to the diagrams shown in Fig. ^ Combined with the 
external quark wave-function renormalization this gives the W^. At tree level, in the massless 
final state limit only is nonzero. Wl can be expressed as 

Wl{yrz) = WliEo,t,z,y)S{z-p)+Wl,iEo,t,z)eiz-p) + WliEo,t,z,XG'), (7) 
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where 

W;„(-B„, t, J, = -16 - ? log' (i^) - J log (i^) log (rt 

+ - ^ + log (i±|) + (2-8^ 

8^. / 2Eot \ 8^. f 2Eot \ 



t "\Eo{t-l) + lJ t " \Eo{l+t) - p, 

1+n , ,„..„^^,„„^l + (^-l)^o^ 



- ^ (log + log (p)) log 



wWf . -2(l + f)(f-p) 2(5i-p) 2(3i + p) 



£;ot2 E^'f^ EoW^ It °\l-t 



Wl{Eo,t,z,XG') = -(l6-^log([±^))log(AG)5(5-p) 



z[z-p + Xg^) h 

+ 0{z-i^p + Xa)f-^^^^±^^^, 
^ ^ t [z-p + Xg') 



and 

/i = ]Jz^ + {XG'-py -2z [Xg' + p) 
h = ^' + 2^0^ (l + f) Xg^+{Xg^-pY -2zp. (8) 

The IR singularities present in the unintegratcd spectra are regulated by a gluon mass Xg in 
this calculation and the divergence cancels between the virtual and bremsstrahlung graphs 
once one integrates over phase space. The divergence directly cancels in integrations of these 
structure functions as the first two bremsstrahlung terms in Wl^ each contribute a factor 
of 81og(AG), while the final term contributes a factor of 8/tlog(^^^ log(AG). For the 
purposes of this paper it is sufficient to numerically integrate the 0{as) spectrum with the 
regulator assigned a small numerical value Xg ~ 10~^. The Li2 {z) functions in Wl are the 
Dilogarithm functions, defined as IA2 (i) = Z^^i or equivalently, as IA2 {z) = /° dt. 
The other structure functions vanish at tree level in the limit mc — > and are IR safe at 
0{as). For these structure functions we find, 
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Wid.z) = Wl(E„,t.z)S(z-p) + Wit(E„.t.z)e(z-p), 
WLiE,„t,z) = -^±ElloJl±l), (9) 



Eot V 1 - t 

W,,[J^O,t,Z) - .^^^^ ^^2^^ ^^3^^ + 

/ 4 2(z + p) (i2_8^_p2) 2z{z + p) z{z'-p') \l 1 + t 
[ Eo' 2Eo' Eo' + 2Eo' J '""^h ~ t> 

W^iy,z) = W^M,t,y)5iz-p) + W^,{Eo,t,z)9iz~p) 



WL[E.,t,y) ^ . ^ log -4 log (p) J (10) 

WUE,,t,z) = -i^ - ^ V/^ + " 7 - 4 ^ ^ 

8 , 4(z-p) (^-p)^-8^ , 5^(^-p) z{z-pf \l^ 

^ I i5^o^ V i^^o' V 2Eo^ Jt^i^Hl-t 

= Wl{E^,t,z)e{z-p) 

H.Wp , - 2(^-p) (^-Pf I 16 (^-P) 2(^-p)^ 

l^4fel-^o,r, - —T'E^A Timi ^ FTs: 



zEaf^ zEoU^ EqH'^ Eo'^t^ 



7{z-p) 3{z-pf 2z{z-p)\l fl + t 



■ V V ^^j^^^Hi^j ^''^ 

W^iy^z) = W^^iE,,t,y)5iz-p) + W,\iE,,t,z)eiz-p) 
^(^(i_p),o,(i±|),4.o,,,)) 

uA(p , ^^ - 11^1^^ 2(lH-3p) , (^-p)^-4p(i-p) , A{z + 3p) 
/8 2(z-p) i-9p i2^2ip-3p2\ 1 /l + t\ 

We have checked the hadron tensor at 0{as) by integrating our results to compare against 
known 0{as) spectra and agree with and the historical [30], but disagree, as do these 
other authors, with .27]. We also agree with the total 0{as) contribution to the decay rate 
in js^. The massless limit of the 0{as) hadron structure functions has been taken for all 



regular terms and we find we agree with the regular terms for a massless final state 



3- 
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C. Lepton Tensor and Phase Space 



To find moments of the triple differential spectrum one must integrate over phase space 
while imposing the experimentally required cut on the lepton energy. With no lepton energy 
cut the phase space is given by the following region [2^, referred to as region Rj, 

l„ /77 ' 1 



^{b-^b^-y) < Ee < -{b+^b^-y) 

(v^+Ag)' < z < (1- V^)' (13) 
< y < 

where b = 1/2 {1 — z + y) and p = ml /ml . Without a lepton energy cut in the phase space, 
the lepton tensor integrated over the lepton energy Eg is 

Lf,u{y, z) =y/3 {-g^^ + q^,qu/y) ■ (14) 

When a minimum cut i?™™ > x mf, is introduced, the phase space and the lepton tensor 
are modified. We do not repeat the derivation of the lepton tensor with a cut here, see 
but note that the phase space splits into three regions when a cut is imposed. These three 
regions correspond to the partitioning of phase space that occurs when the electron energy 
lies below or within the phase space integration range as shown in Fig. 2. We only consider 
cuts below the upper limit of E^ as given in Eq. ()13|). this corresponds to only considering 
cuts where Ef"" < mbil - y^)/2. 

For X < 1/2 (a - ^a-^), where a . 1/2 (l - (^P + Xaf + y) the eut. labelled as 
in Fig. 2, is below the lower limit of the lepton energy in Eq. and as the integration 
over the lepton energy is unaffected, the lepton tensor with this cut reduces to the simple 
expression above with the electron energy integrated over 



-{b-^b^-y)<E,<-{b+^b^-y). (15) 

However, this lepton energy cut still affects the subsequent integrations of z and y by 
imposing the constraint on the range of y , 

2x 

(1 - {^+Xg? - 22^)73^ < (16) 
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z=(l-2i?,)(l-^) 



(Vp + ^g) 




Xi X2 



(\/P + Ag)- 



.•5=(l-2x)(l-^) 




y = (2x)^ 



Rin 
.Rii 



FIG. 2: Phase space diagrams with a lepton energy cut. 



so that the remaining kinematic variables are integrated over the phase space region i?//, 



(17) 



;i-(v^+A 



G) 



2x)- 



2x 



2x 



< y < (1- v^-A 



G) 



For X > 1/2 (a - - v) 

while a; < \/y/2, labelled in the diagram as the cut X2, the phase 
space splits into two regions. The first region Rjjj has the hadron tensor contracted with 
the lepton tensor of Eq. and the range of Ei as given by Eq. ()15j) . The remaining 

phase space variables are then integrated over the range, 



2x)(l-^) < z < [1 



(18) 



{2xy < y < (1-(v/p + Ag) -2x)- 



2x 



2x 



The second region of type Rm combines with the region of phase space where x > 
1/2 (^a — y/a'^ — while x > \/y/2, labelled as xs on Fig. 2. The resulting combined phase 
space has the lepton tensor incorporating the cut within the phase space range and the 
subsequent integration is given by the region Riv, 



X < Ee < -{b+^b^-y), 



(Vp + Xg? < 



< (l-2x)(l-^), 



(19) 



2x 

< y < (l-(v/p + AG)'-2a;) 



l-2x 
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III. HADRONIC MASS MOMENTS 



IS Mass Scheme 



In calculating moments of the hadronic invariant mass spectrum, we use the IS mass and 
the upsilon expansion [2^ . It is well known the pole mass renormalon 



ambiguity 



leads to unnecessarily badly behaved perturbation series for moments when a poor mass 
scheme is chosen. In the IS scheme the renormalon ambiguity is of O (^A.q^qy^ / mfj and so we 
expect the perturbation series of moments of the spectra to be better behaved. We express 
moments in terms of the IS mass, which is related to the h quark pole mass through the 



relation 



22 



m 



IS 



m 



pole 



Ile + — 



11\ o 



4 



(20) 



where £ = log[^ /{rrihasCF)] and Cp = 4/3. The parameter e = 1 determines the order in 
the modified perturbative expansion. Using the upsilon expansion necessitates introducing 
a modified perturbative expansion in order to ensure the cancellation of renormalon ambi- 
guities 2^. When calculating in the IS mass scheme the 0{as^) perturbative corrections 
coming from the mass transformation Eq. (j^Uj) are counted using the parameter e^~^, while 
0(0;^") corrections in the decay rate are counted as e". 

The dependence on the pole mass of the charm quark in our results is eliminated through, 



pole 

— m 



P°i- = mB-rriD + Xi 



2 m/) — mx 
2mT fhr) 



+ Ai 



mx — 2 m^ 
4m|) 



+ AiA 



4m|, — mi, 



is- 



T 



2m|) m\ 



(n + - Pi 



4ml, m^ +^l^QCD 



(21) 



The meson masses mo and ms are the spin averaged meson masses mx = {fnx + 3mx*) /4. 

Aqcd and expand in the parameter Ais 



In this relation we use the fact that 



1 S 

— mt ~ 



A _ "^T IS 

Ai5 = — rrif^ 



(22) 



The perturbative corrections coming from expressing the b quark mass in terms of the 
IS mass are determined by using the definitions of Ais and m^^ and the HQET relationship 
between meson masses and quark masses 

mH = mQ + A — + —r—o 7—0 + ^ ' (23) 



2mc 



4m^ 



4m^ 



m^ 
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where tjih is the hadronic mass, mg is the heavy quark mass, and (i/^- = 3 for the pseudoscaler 
mesons while dn = for the vector mesons. 

In calculating the general hadronic moment, previously calculated moments by Bauer, 
Ligeti, Luke and Manohar (BLLM) were reexamined in the IS mass scheme to check 
results. The results presented in the following sections for the first hadronic moment and 
its variance are different for two reasons. First, in BLLM a l/rfiB expansion was used to 
replace in the expansion of sh in terms of partonic variables. 



4m^ 2m^ 4m^ 4m^ y J 



QCD 



In the results reported in the following sections we always use a 1/mr expansion. The 
corresponding expansion is 

3A2 6A2A1S 3 3 / A^cd\ 

TUB = TUB 5 ^ {T2 + n) H -P2 + C U, . (25) 

rrij; mx mx mx \ rrij. J 

Second, we treat a class of powers of {rfiB — ^)'^ differently than BLLM. When terms are 
generated by changing A = itlb — rrih into A15, 



A = rriB — +Ai5H \ h C e, — 3— , (26) 

2 / mx mx mx V m:^ 



the [rfiB — H|r )" terms are kept only for n < 3. In BLLM this class of terms are treated 
as O (1) although they are formally of order Aq^^j^, leading to this class of higher order 
terms of (m^ — ^)" being kept and summed into the coefficients of the nonperturbative 
parameters. In the following results, the (ffiB — ^)" terms from A are counted as order 
^QCD in the results of Section ITTTl thev are kept only up to 0(^Aqq-q^ in the nonper- 
turbative expansion and up to C(Aqcd) for perturbative terms. Factors of {jfiB — ^) are 
also generated in the replacement of the c and h quark mass and these factors are treated as 
O (1). This implementation of the IS scheme is similar to the the general moment analysis 
of the lepton energy spectrum |2Q|]. The lepton moments are presented in Appendix C and 
can be combined with the hadronic spectrum results to cross check extractions of A15 and 
Ai from these differing spectra. Further hadronic moments that are appropriate for a global 
fit, such as and s%'^ are presented in Appendix B. 
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B. 



Decay Width to ©(a, AqcdM) , ©(A^co/mfe^) 



With the IS scheme implemented as discussed in the previous section, the decay width 
oiB ^X,iu is 



T{B ^ XJu) 



0.5325 - 1.132 — 0.924 — -1.89' 



m-r/2 \m'f/2) \mr/2j 

Ai Xo Ai Ai e Ai q 

- 2.12 ^ - 3.93 + 0.74 - 1.73 \ 

{mr/2y {mr/2y {mr/2y {rnr/2f 

- 5.99 + 4.94 — - 2.98 + 0.99 (27) 

{mr/2f {mr/2f {mr/2f {mr/2f ^ ' 

- 4.96 - 4.94 - 0.080 e + 0.133 e 



+ 0.004 e 



{mr/2Y ' {mr/2Y ' ' mT/2 

Ais V 



^mx/2^ 



where Fq = ^3 ( • Uncertainties in the values of A15, Aj, pi and introduce uncer- 
tainties in the inclusive extractions of |V^fc| using the decay width. In the nonperturbative 
expansion the largest theoretical uncertainty in the extraction of |Vcf,| comes from Ai^ and 
Ai which one can see introduce ~ 2% uncertainties, followed by the higher order nonpertur- 
bative terms which impose an uncertainty of ~ 1% as one can see by examining the results 
for the total decay width and estimating the size of the unknown terms with dimensional 
analysis. 

The size of the OipP^) can be estimated by calculating the /?o contribution to 0{ppy) 



jl?! I33L l34l |. although these terms are not included in this paper. This has been done 
in the IS scheme for a number of observables and the full size of these corrections being 
treated as an error introduces ~ 2% uncertainty in the extraction of |V^b|. The next largest 
uncalculated contributions in the decay width are the 0[as Aj) and 0[as Ai^) terms and the 
^(^QCD /f^b'^c) neglected terms. The size of the 0(^as AQfjjy/mfj terms in the IS scheme 
may be estimated by taking the size of the agAis and multiplying by Aqcd/"^?) ~ 0.1. 
For the first moment, this indicates that the order of the Aj terms is expected to be 
~ 0.01 Aq(jj-, which can be safely neglected in fits to determine the third order parameters in 
the OPE. However, completely uncorrelated uncertainties of this size for both as Aj should be 
used to estimate the error on the fit. The size of the corrections introduced when using the 
mass splitting formula to replace the c quark mass should also be estimated in a fit to extract 
the third order terms in the OPE, as well as uncertainties due to I /ml corrections to the 
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OPE. These number and size of these terms are completely unknown and the uncertainties 
introduced due to these terms can be estimated by introducing completely uncorrelated 
errors of their naive dimensional size. 



C. Integral Hadronic Moments 

The nonperturbative parameters in the decay width can be determined by global fits to 
moments calculated from the decay spectra of B ^ Xciiy such as the hadronic invariant 
mass spectra. The first and second moments of the hadronic invariant mass spectrum have 
been known for some time ^, Is^. The perturbative corrections to these moments were 
obtained by expressing the hadronic moments in terms of the known 0{as) corrections to 
the lepton spectra and the leading order hadronic invariant mass spectra. This technique 
fails when a lepton energy cut is introduced into the phase space, and the general tensor 
results of Section |n] are required. In terms of partonic quantities the hadronic invariant 
mass is defined to be 

sh = {Q - qY = - mBirib {I - z + y) + ml y. (28) 

It is conventional to examine the first hadronic moment once the spin averaged meson 
mass fh\) is subtracted. Moments that give the mean and variance of the hadronic invariant 
mass spectrum with lepton energy cuts of differing values were re-examined recently by 
BLLM. These moments are defined with lepton energy cuts 



nmin 



S, (Er) = {sH - ml) |^^>^.., S, {Er) = {{sH - {sH)f) \E,>Er- ■ (29) 

These moments as functions of the cut on the lepton energy, including the previously 
uncalculated perturbative corrections are as follows. The coefficients stated are for the 
dimensionful nonperturbative parameters listed at the top of the column. The data used in 
the numerical evaluations in this paper are rfiB = 5.3135 GeV, m/j = 1.9730 GeV, rrir = 
9.4603 GeV and the strong coupling is as{mb) = 0.22. For example for the first moment 
with no cut we find 

Si (0) = 0.834 + 1.646 Ais + 0.451 A^^ + 0.16 A^^ + 1.43 Ai - 0.34 As + 0.51 A1A15 
+ 0.07 A2A15 + 0.71 pi - 0.34 p2 + 0.32 n + 0.25 T2 + 0.29 + 0.15 
+ 0.143 e + 0.175 Aise, (30) 
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and in the following tables the leading order term of a moment Si is labelled Sf. 





5? Ais 




A3 


Ai 


A2 


Ai Ais 


A2 Ai5 


Pi 


P2 Ti T2 Ts T4 





0.834 1.646 


0.451 


0.16 


1.43 


-0.34 


0.51 


0.07 


0.71 


-0.34 0.32 0.25 0.29 0.15 


0.5 


0.822 1.623 


0.445 


0.16 


1.44 


-0.30 


0.51 


0.09 


0.72 


-0.34 0.32 0.26 0.29 0.16 


0.7 


0.807 1.592 


0.435 


0.16 


1.46 


-0.24 


0.53 


0.12 


0.75 


-0.34 0.33 0.27 0.29 0.17 


0.9 


0.786 1.549 


0.420 


0.16 


1.51 


-0.14 


0.55 


0.18 


0.79 


-0.34 0.34 0.30 0.30 0.18 


1.1 


0.762 1.496 


0.397 


0.15 


1.57 


0.00 


0.59 


0.26 


0.87 


-0.33 0.35 0.34 0.31 0.21 


1.3 


0.737 1.439 


0.368 


0.15 


1.69 


0.18 


0.66 


0.37 


0.99 


-0.30 0.37 0.41 0.32 0.24 


1.5 


0.719 1.392 


0.334 


0.14 


1.92 


0.42 


0.79 


0.51 


1.23 


-0.23 0.41 0.53 0.33 0.28 



TABLE 1: Coefficients of the nonperturbative parameters for Si (Eq). 





IS 


Contribution 


as Contribution 


Comb 


ined (e) 


77' mill 

Ee 


e 


Ai5 e 


e 


Ai5 e 


e 


Ai5 e 





-0.014 


0.090 


0.157 


0.086 


0.143 


0.175 


0.5 


-0.014 


0.088 


0.143 


0.069 


0.129 


0.157 


0.7 


-0.014 


0.086 


0.139 


0.072 


0.125 


0.159 


0.9 


-0.014 


0.084 


0.134 


0.076 


0.120 


0.160 


1.1 


-0.014 


0.081 


0.128 


0.080 


0.114 


0.161 


1.3 


-0.015 


0.077 


0.121 


0.085 


0.106 


0.162 


1.5 


-0.018 


0.073 


0.117 


0.093 


0.099 


0.166 



TABLE 2: Coefficients of the perturbative parameters for 5*1 {Eq). 



For the S2 (Eq) moments, as explained in Section Qll Al the results differ from those stated 
in BLLM. This difference is formally of higher order, and in the nonperturbative expansion 
the overall effect of the differing implementations of the IS scheme is small. The effect of 
these terms in the perturbative expansion is also small for most moments. However, the 
variance of the hadronic invariant mass spectrum is more sensitive to higher order terms due 
to the cancellation among leading order terms in the nonperturbative expansion. The IS 
scheme as implemented in BLLM found that the variance increased as the lepton energy cut 
was increased due to the dominance of the 0{as) term and the suppression of leading order 
nonperturbative corrections. The 0{as) term is the dominant correction to the variance 
in the rrih 00 limit. The experimentally measured lepton energy cut dependence has 
the variance decreasing as the lepton energy cut increases. When (m_B — ^) terms are 
treated as detailed in Section Fill Al in the IS scheme, the 0{as) term and variance has the 
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experimentally measured dependence on the lepton energy cut, as can be seen in Table 3 
and 4. 





SI 


Ais 






Ai 


A2 


AiAis 


A2 Ais' 


Pi 


P2 


n 


T2 r-s T4 





0.0163 


0.09 


0.08 


-0.03 


-4.87 


1.33 


-1.85 


2.04 


-6.41 


1.51 


-1.04 


-2.87 0.00 0.25 


0.5 


0.0152 


0.08 


0.07 


-0.03 


-4.78 


1.34 


-1.79 


2.07 


-6.48 


1.39 


-1.03 


-2.82 0.00 0.25 


0.7 


0.0153 


0.08 


0.07 


-0.03 


-4.65 


1.35 


-1.70 


2.09 


-6.61 


1.23 


-1.01 


-2.77 0.00 0.25 


0.9 


0.0167 


0.09 


0.07 


-0.04 


-4.48 


1.37 


-1.56 


2.12 


-6.84 


0.99 


-0.98 


-2.69 0.00 0.25 


1.1 


0.0195 


0.10 


0.08 


-0.04 


-4.26 


1.40 


-1.36 


2.17 


-7.23 


0.68 


-0.95 


-2.61 0.00 0.24 


1.3 


0.0224 


0.11 


0.09 


-0.04 


-4.01 


1.44 


-1.10 


2.25 


-7.88 


0.28 


-0.92 


-2.53 0.00 0.24 


1.5 


0.0227 


0.11 


0.09 


-0.04 


-3.76 


1.52 


-0.77 


2.38 


-9.05 


-0.22 


-0.91 


-2.48 0.00 0.24 



TABLE 3: Coefficients of the nonperturbative parameters for 5*2 {Eq 





IS 


Contribution 


Og Contribution 


Combined (e) 


771 mill 

Ee 


e 


Ais e 


e 


Ai5 e 


e 


Ai5 e 





0.105 


0.163 


0.551 


0.424 


0.656 


0.588 


0.5 


0.102 


0.159 


0.163 


-0.202 


0.265 


-0.043 


0.7 


0.099 


0.154 


0.099 


-0.266 


0.198 


-0.112 


0.9 


0.094 


0.147 


0.057 


-0.283 


0.151 


-0.136 


1.1 


0.089 


0.140 


0.029 


-0.273 


0.118 


-0.132 


1.3 


0.085 


0.135 


0.010 


-0.240 


0.096 


-0.112 


1.5 


0.083 


0.133 


0.002 


-0.211 


0.083 


-0.078 



TABLE 4: Coefficients of the perturbative parameters for S2 (Eq). 

The moments Si and 5*2 with a lepton energy cut of 1.5 GeV have been experimentally 
measured by CLEO l3 |. 



^1 (1.5 GeV) = 0.251 ± 0.066 GeV^ 
^2 (1.5 GeV) = 0.576 ± 0.170 GeV^ 



(31) 
(32) 



Using this data, we can extract values of A15 and Ai, for comparison with extractions from 
the lepton energy spectrum. We estimate the theoretical error on the extraction due to the 
unknown third order terms in the usual way [91 using the results of recent ffis when they 
substantially improve our knowledge of these terms beyond dimensional analysis. We use 
the HQET vector pseudoscalar mass splitting constraint to determine A2 = 0.12GeV^ and 
the mass splitting formula to third order 
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P2-T2- Ti 



K, (mc) ml iric Am b — mb m? Am o 



(33) 



rub — rricK [mc) 

to reduce the number of free parameters. A positive value of pi is imposed in accordance 
with vacuum saturation 3^ and is drawn from the range [0,0.125] GeV^. The unknown 
matrix elements are then drawn from a flat distribution, the unknown third order terms 
are drawn from between ±0.125 GeV^ while Ai in drawn from [—600, 0] MeV^ in accordance 
with its full constrained range from the BLLM fit. We then extract the following values for 
Ai and A15: 



Ai5=[_0.13 ± 0.05, ± 0.09„3] GeV 
Ai = [-0.24 ± 0.02, ± 0.09„3] GeV^ 



(34) 



The perturbative errors are estimated by using the two loop running of as to vary the scale 
of as (/i) between m;,/2 < fi < 2mb- Adding the theoretical errors in quadrature we obtain 



m, 



IS 



4.86 ± 0.10 GeV, which is in excellent agreement with the inclusive extraction using 



the lepton energy moments, m]^ = 4.84 ±0.10 GeV [20| and in agreement with the results of 
BLLM where m^^ = 4.74 ±0.10 GeV, despite the differences in the expansion and the larger 
number of observables in the fit. This extracted IS mass m^ = 4.86 ± 0.10 GeV translates 
into a value of the MS mass ffihimh) = 4.34 ± 0.09 GeV which can be compared with other 



extractions of the MS mass 
bb production cross section 



4.25 ±0.08 GeV 



4l|. 



36 



21 1 such as the MS mass found by examining moments of the 



37, 



39l I40I 1 a recent analysis of which finds mb{mb) 



D. Fractional Moments 

1. The l/rric Expansion 



In integer moments such as Si and 5*2, a Aqcd/'o^c expansion only enters the predicted 
value of a moment through the mass transformation relationship Eq. ( |^ ). Fractional 
moments have additional cuts in the complex q ■ v plane due to the branch cut starting at 
sh = when sh is taken to a noninteger power. These branch cuts are separated from the 



physical cut by a scale set by mc, as the physical branch cut begins at sh 



m 



D- 



As mr 







these cuts coalesce and one would expect predictions for fractional moments in this limit to 
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be ill defined, as discussed in recent work [19]. We find an explicit mhIs.QCD/m1 expansion 
in calculations of fractional moments of sh'-, the neglected terms in this expansion are nu- 
merically suppressed for hadronic invariant mass observables leading to a small numerical 
uncertainty being introduced. This can be shown by examining a general moment of the 
squared hadronic invariant mass s^. The dependence of the general moment as a function 
of n is found by performing an expansion of s^, 

The coefficient functions Cf are O {j\^Qj^/ml^. Expanding up to O {jVq^Qj^/m^ in the 
nonperturbative expansion we find, 



^^"^ [Cl + J {Cl + yd) + {d + yCl + f d 

+ (pl + yd + f Cl + f Cl) ] , (36) 

where the are functions of n and the nonperturbative matrix elements ^. For integer 
moments this expression has no 1/z dependence. However, for non-integer moments in this 
range one obtains contributions of order where A; > n is the ceiling of the fractional 
moment power n. As the lower limit of 2; is p = ml/ml, this corresponds to a rribAQCD/^c 
expansion entering into the calculations of fractional moments. This expansion does not 
seem to introduce a large uncertainty for fractional moments compared to integer moments 
as the neglected class of terms are numerically suppressed for n values in the range [0, 3] 
but the best way to estimate the uncertainty introduced is under study. In the following 
investigation of hadronic fractional moments no additional uncertainty is added to account 
for this theoretical error and we examine how known sources of error primarily from unknown 
matrix elements can be reduced. 

When examining a general moment to obtain interesting observables, we expand in the 
ratio Aqcd/''^q and then examine the n and i?™™ dependence of the coefficient functions of 
the nonperturbative matrix elements. The essential observation motivating this approach is 
that one is allowed to choose n and Ef^^'^ within a range of reasonably accessible experimental 



^ These C* coefScient functions are reported in the Appendix. 
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values, in order to maximize the utility of a measured moment in obtaining information on 
the nonperturbative matrix elements. We define a general moment function, 

S[n,E,„m,Ee,]= \'l\ (37) 

\Ei>E2 



so that 



nminl —2 



Si (^Ef'') = S[l, ^ 0, E^ 

[Ef'') = S% Ef", 0, Ef''] - {S[l, 0, E, 



2 



Our search of the hadronic mass moments is restricted to the parameter space, 

m < 3 , n < 3 , 0.5 GeV < E^^f < 1.5 GeV, (38) 

to ensure a well behaved OPE. In this parameter space we find two types of moments of 
interest, moments that allow the OPE to be precisely tested for deviations from experiment 
and moments that allow one to extract the IS mass with minimal error. We consider each 
type of moment in the following sections. 

2. OPE Testing Moments 

A discrepancy of the prediction of the OPE when compared with data can come from a 
number of possible sources when one is considering percent level extractions of \ Vcb\'- a higher 
order matrix element that is being neglected could be anomalously large, the OPE itse 



45 



42 



f could 



not be converging or quark-hadron duality violation could effect determinations 

finding moments where only the leading order unknown nonperturbative 



parameters are suppressed and checking the predicted values against experiment, one can 
assess the theoretical error that is being assigned to the inclusive extraction of \ Vcb\ in a clear 
and unambiguous fashion. This technique ||2|| has recently been used to test the OPE in the 



lepton energy spectrum. Measurements of these OPE testing observables in this spectrum 
indicated that the OPE is a valid description of the data to the percent level j4S|. It is 
important to note that the OPE testing moments presented allow one to check that the 
error assigned in the extractions of \Vcb\ and m^, is large enough to account for all of these 
possible effects if the OPE testing moments properly describe the data 
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A selection of these OPE testing moments for the hadronic invariant mass spectrum is 
presented in Table 5. The nonperturbative parameters for these moments are not suppressed 
except for the leading unknown nonperturbative terms A15 and Ai. A typical OPE testing 
moment is as follows, 



Dl = 5[2, 0.5, 2.2, 0.7] 
= 0.7779 



Ai5 / Ais V ( Ais \^ Ai 

1 + 0.168 — + 0.345 — ^- +0.501 — ^- +0.01 



mr/2 ■ \mr/2j ' \mr/2 J ' {mr/2y 



Ao Ai Ai q Ao Ai Q jOi 

-1.05- ^-1.94 ^ ^% -8.75 ^ ^^3 +6.18^ 



{mr/2y {mr/2y {mr/2y {mr/2y 



-0.25 — ^^ + 0.27 ^^ + 3.71 ^^-1.50- 



{mr/2y {mr/2y {mr/2y {mr/2y 



-1.61 - 0.0094 e - 0.008 e 

{mr/2y mr/2 



(39) 



Varying the unknown parameters in the same way, and treating the leading order nonper- 
turbative parameters as unknowns and varying them over the region Ai^ = —0.13 + 0.1 GeV, 
Ai = —0.3 ± 0.3 GeV^, this moment is predicted to be 



Dl = 0.7686 + 0.0018e + 0.0040n.p. . (40) 

The perturbative error is obtained by the scale variation in the standard way. As in the 
lepton spectra, the OPE testing moments are such that with no nonperturbative input other 
than the known value of A2, wc can predict the value of a moment to an accuracy of 1%. 

As the error on the nonperturbative terms is reduced with global fits, it is important to 
cross check with the predicted and measured values of these moments in order to ensure the 
error on \Vcb\ and is not being underestimated. Table 5 and 6 present a selection of 
OPE testing moments in the parameter space examined. 







Ais 




A3 


Ai 


A2 


AiAis 


A2 Ai5 


Pi 


P2 


n 


T2 


T3 




Dl 


0.7779 


-0.028 


0.01 


0.00 


0.00 


-0.04 


-0.01 


-0.06 


0.05 


-0.00 


0.00 


0.03 


-0.01 


-0.01 


Dl 


0.8845 


-0.019 


0.00 


0.00 


0.00 


-0.02 


-0.01 


-0.03 


0.03 


0.00 


0.00 


0.01 


-0.01 


-0.01 


Dl 


0.7829 


0.047 


0.10 


0.04 


0.00 


-0.13 


-0.03 


-0.16 


0.06 


0.03 


0.00 


0.03 


-0.02 


-0.02 


Dj 


1.9030 


0.166 


-0.08 


-0.03 


0.00 


0.22 


0.08 


0.38 


-0.27 


0.00 


-0.01 


-0.15 


0.06 


0.07 


Dl 


2.428 


0.040 


-0.46 


-0.20 


0.00 


0.66 


0.19 


0.99 


-0.62 


-0.08 


-0.03 


-0.29 


0.12 


0.14 
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TABLE 5: Coefficients of the nonperturbative parameters for OPE testing Moments. The 
OPE testing moments in the table are defined in the following way, D\ — S [2.0, 0.5, 2.2, 0.7], 
Dl = 5 [1.9, 0.6, 2, 0.7] , Dl = 5 [2.6, 0.6, 2.9, 1], Dj = 5 [2.4, 1, 1.9, 0.8] and = 
5[2.9, 1.4,2.2, 1.3]. 





IS as'^ Term 


as Term 


Combined O (e) 


Predicted Value for Moment 


Label 


e Ai5 e 


e Ai5 e 


e Ai5 e 




Dl 
Dl 
Dl 
Dj 

d: 


-0.001 -0.005 
-0.001 -0.003 
-0.001 -0.004 
0.006 0.029 

0.010 0.050 


-0.006 -0.001 
-0.005 -0.004 
-0.005 -0.012 
0.025 -0.012 

0.026 -0.036 


-0.007 -0.007 
-0.006 -0.007 
-0.006 -0.016 
0.030 0.017 

0.036 0.011 


0.7686 ±0.0018 e ± 0.0040(N.P.) 
0.8804 ±0.0014 e ± 0.0025(N.P.) 
0.7582 ±0.0016 e ± 0.0067(N.P.) 
1.9448 ±0.0078 e ± 0.026(N.P.) 

2.5380 ±0.0102 f ±0.017(N.P.) 



TABLE 6: Coefficients of the perturbative parameters for OPE testing moments and their 
predicted value. 



3. Moments to measure ml^ with minimal error 



Moments that allow a direct measurement of the b quark mass with minimal dependence 
on the unknown Ai nonperturbative parameter have also been found. Moments of this 
type are important as a measurement of m^^ to this precision will be an important step in 

reducing the error on \ Vuh\ as well as extracting \ Vcb\- These moments are particularly suited 
to being used in a fit to extract the IS mass. An example of this type of moment is 



Bl = 5[3, 0.5, 0.5, 0.9] 
= 51.9318 



Ais / Ais- 

1± 3.215 — i^± 2.609 ' * 



mT/2 



± 3.67 



A, 



0.004 



{mr/2y 



±6.14 



AiA 



IS 



imr/2y 



-29.15 



imr/2y 

^2 



± 45.49 



mx/2 J 
A2 Ais 



2.216 



A 



15 



{mr/2y 



+ 13.26 



imr/2y 



27.78 



mx/2 J 
Pi 



±0.40 



Ai 



(mT/2)' 



± 10.64 



(mT/2)' 



±4.10 



(mT/2)^ 

P2 



{mr/2y 



(mT/2)- 



± 0.084 e± 0.154 e 



A 



IS 



my/2 



(41) 



Note that moments of this type have a strong dependence on A15 and a weak dependence 
on Ai while the coefficients of the higher order terms in the nonperturbative series are not 
suppressed. These moments are thus well suited to measure the b quark mass with minimal 
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error as the largest source of theoretical error is suppressed in a controlled fashion. Esti- 
mating the error on the extraction in the usual way one finds the mass extracted from 
this moment will have a theoretical error due to unknown matrix elements and perturba- 
tive terms of ±50MeV(N.P.) ± 3MeV(e), where the error is dominated by the unknown 
nonperturbative corrections at third order. Adding these errors in quadrature one obtains 
a theoretical error in the extraction ~ 50 MeV. It is important to experimentally measure 
moments of this type for a precise value of the b quark mass to be extracted from this 
spectrum. This error assessment in the extraction of is assigned in accordance with 
how theoretical error is assigned in the OPE testing moments. By measuring both the OPE 
testing moments and the moments presented in this section, one can extract m^^ with an ex- 
perimentally tested theoretical error; comparisons of the b quark mass extracted in this way 
with extractions using other techniques will be a useful cross check of theoretical techniques. 
A selection of moments of this type is given in Table 7 and 8. 







Ais 






Ai 


A2 


AiAis 


A2 Ais 


Pi 


P2 


Tl 


T2 Ta T4 


Bl 


51.9318 35.300 


6.06 


-1.09 


0.93 


8.51 


3.01 


22.32 


-13.63 


2.01 


0.00 


-14.30 6.51 5.22 


Bl 


10.2470 


4.850 


1.38 


0.37 


0.02 


1.74 


0.43 


3.25 


-2.84 


-0.03 


-0.04 


-2.03 0.88 0.72 


Bl 


12.2674 


6.116 


1.48 


0.28 


0.31 


1.35 


0.79 


3.62 


-3.50 


0.24 


0.00 


-2.43 1.10 0.88 


B\ 


4.7852 


1.659 


0.55 


0.20 


0.10 


0.55 


0.16 


0.93 


-0.95 


-0.11 


0.01 


-0.62 0.30 0.23 



TABLE 7 : The nonperturbative parameters of moments to measure the IS mass accurately. 
In the table above the moments are defined in the following way B\ — [3, 0.5, 0.5, 0.9], 
Bl = S [2.4, 0.5, 1.2, 1.3], Bl = S [2.5, 0.5, 1, 1] and Bj = S [2, 0.5, 1.3, 1.3]. 





IS tts^ Term 


as Term 


Combined (e) 


Error in extraction of IS mass 




e Ais e 


e Ais e 


e Ai5 e 




Bl 
Bl 
Bl 
B', 


-0.003 2.231 
0.042 0.353 
0.042 0.454 
0.017 0.121 


4.345 -0.536 
0.383 -0.336 

0.564 -0.318 
0.090 -0.158 


4.343 1.70 
0.426 0.017 
0.606 0.137 
0.107 -0.004 


±3MeVe ±50MeV(N.P.) 
±5MeVe ±55MeV(N.P.) 
ilMcVe ±54MeV(N.P.) 
±2MeVe ±59MeV(N.P.) 



TABLE 8: The perturbative parameters of moments to measure the IS mass accurately, 
and the error estimated in using the moment to extract the IS mass. 
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IV. CONCLUSIONS 



We have presented the 0{as) corrections to the structure functions of the hadronic tensor 
for S — > Xc^P . The 0{as) and ©(cts Aqcd/"t.6) perturbative corrections for lepton energy 
moments and hadronic invariant mass moments have been calculated. The effects of the 
charm quark expansion in fractional moments was shown to be small and moments that allow 
one to extract the nonpcrturbative parameters rcvelant for a percent level determination of 
\Vcb\ from inclusive B decay were presented. Using the techniques outlined, a b quark 
mass measurement with a theoretical error at the 50MeV or better should be possible 
using the inclusive semileptonic decay data. This theoretical error assessment, including 
the assumption of negligible quark-hadron duality violation that this analysis relies upon, 
is directly testable with the OPE testing moments presented in this paper. Fits based on 
the results presented should allow an extraction of |V"cfe| with ~ 2% theoretical error. As 
the lepton energy cut dependence of the O {ag Aqcb/ ttii,) term is now known, the largest 
estimated theoretical uncertainty in inclusive extractions of \Vcb\ and ml^ comes from the 
0{al) corrections. 
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APPENDIX A: sh EXPANSION 

The coefficient functions of the general hadronic moment sh in terms of the pole mass 
and the HQET local operator expansion are as follows: 



o_i riA n (Ai + 3 As - (n - 1)A2) (n - 1) n A ((n - 2) - 3 (Al + 3 A2)) 
nib 2 m^^ D rrib'^ 



n (ti + 3 T2 + Ta + 3 r4 - pi - 3 ps) 



3 



(Al) 



i_ A (Al +3A2 -2nA2) nA ((n-1) A2-2 (A1 + 3A2)) 

^0 — + n 2 ^ o 3 

rrib 2 nib 2 nib 

(ti + 3 T2 + T3 + 3 r4 - pi - 3 P2) 



Co = + ^ + ' , (A3) 
2 nib 2 nib^ 



A nib 

A (nA2 - Al -3A2) 
2 nib'^ 2 nib 

A3 

Snib^ 

A_ _ (2 A^ + Ai + 3A2) _ (n - 1) A {{n - 2) A^ - 2 (Ai + 3 A2)) 

nib 2 nib"^ 2 m^^ 

(ti + 3 r2 + T3 + 3 r 4 - pi - 3 P2) 

Anib^ 

A2 2A ((n-1) A2-A1-3A2) 



(A2) 



_^ V'l-ru/2-r'3^^'4 — ^1 — ^^2; ^^g^j 



Ci^ = '-^ ^ 3 (A6) 

A^ 

A^^A((„-2)A^ A.-3A.) 
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APPENDIX B: HADRONIC INVARIANT MASS MOMENTS FOR FIT 



pmin 




A2 

^HS 


A?5 


Ai 


A2 


AiAi5 


A2 Ai5 


Pi 


P2 n T2 T3 r4 





2.1799 0.411 


0.133 


0.07 


0.39 


-0.05 


0.07 


0.06 


0.35 


-0.11 0.10 0.11 0.07 0.03 


0.5 


2.1771 0.406 


0.131 


0.07 


0.39 


-0.04 


0.07 


0.07 


0.35 


-0.11 0.10 0.11 0.07 0.03 


0.7 


2.1735 0.399 


0.129 


0.07 


0.40 


-0.03 


0.07 


0.07 


0.35 


-0.11 0.10 0.11 0.07 0.04 


0.9 


2.1685 0.388 


0.125 


0.07 


0.40 


-0.01 


0.08 


0.08 


0.37 


-0.10 0.10 0.12 0.08 0.04 


1.1 


2.1625 0.375 


0.119 


0.06 


0.42 


0.02 


0.09 


0.10 


0.39 


-0.09 0.10 0.13 0.08 0.05 


1.3 


2.1565 0.362 


0.112 


0.06 


0.44 


0.07 


0.10 


0.12 


0.43 


-0.08 0.11 0.14 0.08 0.06 


1.5 


2.1522 0.350 


0.104 


0.06 


0.50 


0.12 


0.12 


0.14 


0.51 


-0.05 0.12 0.17 0.08 0.07 



TABLE 9: Coefficients of nonperturbative parameters for 5*1/2 = >S'[0.5, E^/^^^, 0, E'™"']. 





IS 


Contribution 


as Contribution 


Combined (e) 


^min 


e 


Ai5 e 


e 


Ai5 e 


e 


Ai5 e 





-0.003 


0.023 


0.026 


0.006 


0.023 


0.029 


0.5 


-0.003 


0.022 


0.031 


0.014 


0.027 


0.037 


0.7 


-0.003 


0.022 


0.031 


0.017 


0.028 


0.039 


0.9 


-0.003 


0.021 


0.031 


0.020 


0.027 


0.041 


1.1 


-0.003 


0.020 


0.030 


0.021 


0.027 


0.041 


1.3 


-0.004 


0.019 


0.029 


0.023 


0.025 


0.042 


1.5 


-0.004 


0.018 


0.028 


0.025 


0.024 


0.043 



TABLE 10: Coefficients of perturbative parameters for ^1/2 = ^[0.5, Ef"", 0, Ef""]. 





Ai5 




A?5 


Ai 


A2 


Ai Ai5 


A2 Ai5 


Pi 


P2 n T2 T3 T4 





10.286 5.377 


1.833 


0.56 


3.72 


-1.00 


1.88 


0.21 


0.63 


-0.72 0.76 0.21 0.85 0.50 


0.5 


10.249 5.296 


1.801 


0.55 


3.77 


-0.87 


1.90 


0.30 


0.67 


-0.74 0.77 0.24 0.86 0.51 


0.7 


10.200 5.189 


1.757 


0.54 


3.85 


-0.67 


1.95 


0.42 


0.72 


-0.77 0.78 0.29 0.87 0.55 


0.9 


10.133 5.040 


1.689 


0.53 


4.00 


-0.35 


2.04 


0.63 


0.81 


-0.81 0.81 0.37 0.88 0.60 


1.1 


10.054 4.859 


1.597 


0.51 


4.25 


0.09 


2.19 


0.92 


0.97 


-0.84 0.85 0.50 0.91 0.67 


1.3 


9.976 4.669 


1.485 


0.48 


4.66 


0.68 


2.46 


1.32 


1.22 


-0.84 0.92 0.72 0.94 0.76 


1.5 


9.919 4.515 


1.366 


0.44 


5.40 


1.43 


2.96 


1.83 


1.69 


-0.72 1.05 1.10 0.98 0.89 



TABLE 11: Coefficients of nonperturbative parameters for ^3/2 = ^[1.5, Ef"", 0, Ef""]. 
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IS 


Contribution 


as Contribution 


Comb 


ined (e) 




e 


Ms e 


e 


Ms e 


e 


Ms e 





-0.041 


0.264 


0.603 


0.402 


0.562 


0.667 


0.5 


-0.041 


0.261 


0.483 


0.225 


0.443 


0.485 


0.7 


-0.041 


0.255 


0.458 


0.215 


0.417 


0.470 


0.9 


-0.041 


0.248 


0.432 


0.216 


0.391 


0.464 


1.1 


-0.042 


0.239 


0.405 


0.222 


0.363 


0.462 


1.3 


-0.045 


0.229 


0.380 


0.236 


0.335 


0.465 


1.5 


-0.053 


0.215 


0.363 


0.263 


0.310 


0.478 



TABLE 12: Coefficients of perturbative parameters for ^3/2 = S[1.5, Ef"", 0, Ef% 





sL Ms 




A3 


Ai 


A2 


AiAis 


A2 Ais 


Pi 


p2 n 


r2 n n 





22.297 15.253 


5.984 


1.66 


7.93 


-1.88 


5.00 


1.62 


-0.89 


-1.12 1.47 


-0.89 2.23 1.41 


0.5 


22.189 15.00 


5.856 


1.64 


8.11 


-1.50 


5.10 


1.89 


-0.86 


-1.25 1.50 


-0.79 2.25 1.47 


0.7 


22.048 14.675 


5.685 


1.60 


8.39 


-0.94 


5.28 


2.28 


-0.80 


-1.41 1.55 


-0.63 2.28 1.55 


0.9 


21.857 14.222 


5.438 


1.54 


8.88 


-0.06 


5.60 


2.91 


-0.69 


-1.63 1.63 


-0.38 2.32 1.68 


1.1 


21.634 13.684 


5.120 


1.46 


9.65 


1.18 


6.13 


3.82 


-0.50 


-1.87 1.76 


0.01 2.38 1.87 


1.3 


21.413 13.122 


4.749 


1.36 


10.90 


2.82 


7.04 


5.03 


-0.16 


-2.05 1.97 


0.63 2.46 2.12 


1.5 


21.252 12.671 


4.378 


1.24 


13.07 4.93 


8.69 


6.61 


0.48 


-2.00 2.30 


1.64 2.57 2.44 



TABLE 13: Coefficients of nonperturbative parameters for S2a = 5'[2, E'™™, 0, E'™™]. 





IS 


Contribution 


Os Contribution 


Combined (e) 




e 


Ms e 


e 


Ms e 


e 


Ms e 





-0.104 


0.689 


1.929 


1.367 


1.826 


2.045 


0.5 


-0.103 


0.679 


1.425 


0.596 


1.321 


1.275 


0.7 


-0.103 


0.666 


1.322 


0.536 


1.220 


1.202 


0.9 


-0.104 


0.647 


1.227 


0.526 


1.124 


1.173 


1.1 


-0.107 


0.624 


1.135 


0.541 


1.028 


1.164 


1.3 


-0.115 


0.596 


1.051 


0.578 


0.937 


1.174 


1.5 


-0.135 


0.561 


0.995 


0.656 


0.860 


1.217 



TABLE 14: Coefficients of perturbative parameters for S2a (Eq). 



APPENDIX C: LEPTON ENERGY MOMENTS FOR FIT 

Lepton energy moments appropriate for extracting Ais and Al from previous work [2^ 
in terms of the inverse upsilon mass expansion and A15 = ^ — definition of Eq. ()22|) 
via A'fl'*"'^ = (rfiB — + Ais are as follows. The general moment is defined for the lepton 
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spectrum in an identical fashion to the general moment for the hadronic invariant mass 
spectrum. 



R[n, Ee, , m, E,,] - , (CI) 



^min 




Al5 






Ai 


A2 


AiAi5 


A2 


Pi 


P2 




T"2 


^3 







1.3920 


-0.075 


-0.02 


0.00 


-0.10 


-0.21 


-0.04 


-0.05 


-0.03 


0.01 


-0.02 


-0.01 


-0.03 


-0.03 


0.5 


1.4216 


-0.074 


-0.02 


0.00 


-0.10 


-0.21 


-0.04 


-0.05 


-0.03 


0.01 


-0.02 


-0.01 


-0.03 


-0.03 


0.7 


1.4611 


-0.073 


-0.02 


0.00 


-0.10 


-0.20 


-0.04 


-0.05 


-0.03 


0.01 


-0.02 


-0.02 


-0.03 


-0.03 


0.9 


1.5173 


-0.073 


-0.02 


0.00 


-0.10 


-0.20 


-0.04 


-0.05 


-0.04 


0.01 


-0.02 


-0.02 


-0.03 


-0.03 


1.1 


1.5884 


-0.073 


-0.02 


0.00 


-0.10 


-0.19 


-0.04 


-0.05 


-0.04 


0.00 


-0.02 


-0.02 


-0.02 


-0.03 


1.3 


1.6724 


-0.074 


-0.02 


0.00 


-0.10 


-0.18 


-0.04 


-0.05 


-0.04 


0.00 


-0.02 


-0.02 


-0.02 


-0.03 


1.5 


1.7674 


-0.076 


-0.02 


0.00 


-0.11 


-0.17 


-0.04 


-0.05 


-0.05 


-0.01 


-0.02 


-0.03 


-0.02 


-0.03 



TABLE 15: Coefficients of the nonperturbative parameters for Ri = i?[l, Ef"", 0, Ef% 





IS 






as 


a 


s^(3o 


Combined Terms 


TTTnin 
^/ 


e 


Ai5 e 




e 


Ai5 e 




e 


Ais e 







0.004 


0.001 


0.006 


-0.001 


-0.001 


-0.002 


0.003 


0.000 


0.003 


0.5 


0.004 


0.001 


0.006 


-0.001 


-0.001 


-0.003 


0.003 


0.000 


0.002 


0.7 


0.004 


0.001 


0.006 


-0.001 


-0.001 


-0.004 


0.002 


0.000 


0.002 


0.9 


0.004 


0.001 


0.006 


-0.002 


-0.001 


- 0.005 


0.002 


0.000 


0.001 


1.1 


0.004 


0.001 


0.006 


-0.002 


-0.001 


-0.006 


0.001 


0.000 


0.000 


1.3 


0.004 


0.001 


0.006 


-0.003 


-0.001 


-0.007 


0.001 


0.000 


-0.001 


1.5 


0.004 


0.001 


0.006 


-0.003 


-0.001 


-0.007 


0.001 


0.000 


-0.001 



TABLE 16: Coefficients of the perturbative parameters for Ri {Eq). 



Ee 


vr 


Ai5 




A?,s 


Ai 


A2 


Ai Ais 


A2 Ais 


pi 


P2 


Tl 


r2 


^3 


^4 





0.1804 


-0.032 


0.00 


0.00 


-0.05 


-0.07 


-0.01 


-0.01 


-0.04 - 


0.01 


-0.01 


-0.02 


-0.01 


-0.01 


0.5 


0.1542 


-0.032 


0.00 


0.00 


-0.05 


-0.06 


-0.01 


-0.01 


-0.04 - 


0.01 


-0.01 


-0.02 


-0.01 


-0.01 


0.7 


0.1280 


-0.030 


0.00 


0.00 


-0.05 


-0.06 


-0.01 


-0.01 


-0.04 - 


0.01 


-0.01 


-0.02 


-0.01 


-0.01 


0.9 


0.0988 


-0.028 


0.00 


0.00 


-0.05 


-0.06 


-0.01 


-0.01 


-0.04 - 


0.01 


-0.01 


-0.02 


-0.01 


-0.01 


1.1 


0.0705 


-0.025 


0.00 


0.00 


-0.05 


-0.05 


-0.01 


-0.01 


-0.04 - 


0.01 


-0.01 


-0.02 


-0.01 


-0.01 


1.3 


0.0458 


-0.021 


0.00 


0.00 


-0.05 


-0.04 


-0.01 


-0.01 


-0.04 - 


0.01 


-0.01 


-0.02 


-0.01 


-0.01 


1.5 


0.0261 


-0.017 0.00 


0.00 


-0.04 


-0.03 


-0.01 


-0.02 


-0.04 - 


0.01 


-0.01 


-0.02 


-0.01 


-0.01 



Table 17: Coefficients of the nonperturbative parameters for variance 

Vi = {R[2, Ef"^, 0, Ef''] - R[l, Ef'\ 0, Ef'^f). 
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IS 


Q 




Ois 






Combined Terms 




e 


Ai5 e 


^BLM 


e 


Ai5 e 


^BLM 


e 


Ai5 e 


^BLM 





0.002 


0.000 


0.003 


-0.002 


0.000 


-0.002 


0.000 


0.000 


0.001 


0.5 


0.002 


0.000 


0.003 


-0.002 


-0.001 


-0.001 


0.000 


0.001 


0.001 


0.7 


0.002 


0.000 


0.003 


-0.002 


0.000 


-0.001 


0.000 


0.000 


0.002 


0.9 


0.001 


0.000 


0.002 


-0.001 


0.000 


- 0.001 


0.000 


0.000 


0.002 


1.1 


0.001 


0.000 


0.002 


-0.001 


0.000 


-0.001 


0.000 


0.000 


0.002 


1.3 


0.001 


0.000 


0.002 


-0.001 


0.000 


-0.001 


0.000 


0.000 


0.001 


1.5 


0.001 


0.000 


0.002 


-0.001 


0.000 


0.000 


0.000 


0.000 


0.001 



TABLE 18: Coefficients for lepton variance perturbative parameters. 



pmin 




Ai5 




A?, 


Ai 


A2 


Ai Ai5 


A2 Ai5 


Pi 


P2 




T2 


T3 


T4 





-0.0376 


0.001 


0.002 


0.0 


-0.01 


0.02 


0.0 


0.01 


-0.03 


-0.01 


0.0 


-0.01 


0.0 


0.0 


0.5 


-0.0207 


0.0 


0.002 


0.0 


-0.01 


0.01 


0.0 


0.0 


-0.03 


-0.01 


0.0 


-0.01 


0.0 


0.0 


0.7 


-0.0105 


0.0 


0.001 


0.0 


-0.01 


0.01 


0.0 


0.0 


-0.03 


-0.01 


0.0 


-0.01 


0.0 


0.0 


0.9 


-0.0036 


-0.002 


0.001 


0.0 


-0.01 


0.01 


0.0 


0.0 


-0.03 


-0.01 


0.0 


-0.01 


0.0 


0.0 


1.1 


-0.0001 


-0.002 


0.001 


0.0 


-0.01 


0.0 


0.0 


0.0 


-0.02 


-0.01 


0.0 


-0.01 


0.0 


0.0 


1.3 


0.0009 


-0.002 


0.0 


0.0 


-0.01 


0.0 


0.0 


0.0 


-0.02 


-0.01 


0.0 


-0.01 


0.0 


0.0 


1.5 


0.0009 


-0.001 


0.0 


0.0 


-0.01 


0.0 


0.0 


0.0 


-0.02 


0.0 


0.0 


0.0 


0.0 


0.0 



TABLE 19: Coefficients for the nonperturbative parameters of 

V2 =< {R[l,Ef'\0,Ef'']- < R[l,Ef'\0,Ef''] >f >. 





IS Q 




s^/5o 


as 




Combined Contributions 


^min 


e 


Ai5 e 


^BLM 


e Ai5 e 


^BLM 


e 


Ai5 e 


^'blm 





0.000 


0.000 


0.000 


0.000 0.001 


0.000 


0.000 


0.001 


0.000 


0.5 


0.000 


0.000 


0.000 


0.001 0.000 


0.000 


0.001 


0.000 


0.000 


0.7 


0.000 


0.000 


0.000 


0.000 0.000 


0.000 


0.000 


0.000 


0.000 


0.9 


-0.008 


0.000 


0.000 


0.007 0.001 


0.000 


-0.001 


0.001 


0.000 


1.1 


-0.007 


0.000 


0.000 


0.005 0.000 


0.000 


-0.002 


-0.001 


0.000 


1.3 


-0.006 


0.000 


0.000 


0.005 0.000 


0.000 


-0.001 


0.001 


0.000 


1.5 


-0.005 


0.000 


0.000 


0.003 0.000 


0.000 


-0.002 


0.000 


0.000 



TABLE 20: Coefficients for the perturbative parameters of V2 
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Al5 




A?,, 


Ai 


A2 


A1A15 


A2 Ai5 


Pi 


P2 




T2 


T3 T4 




0.5452 


0.001 


-0.003 


-0.01 


0.002 


-0.02 


0.00 


-0.01 


0.01 


0.01 


0.00 


0.01 


0.01 0.00 


D2 


1.7626 


0.014 


0.014 


0.00 


0.001 


0.09 


0.01 


0.05 


-0.01 


-0.02 0.00 


-0.01 0.01 0.01 


D3 


0.5215 


-0.011 


-0.009 


0.00 


-0.002 


-0.04 


-0.01 


-0.03 


0.01 


0.01 


0.00 


0.00 


0.00 0.00 




0.6051 


-0.015 


-0.011 


-0.003 


-0.006 


-0.04 


-0.01 


-0.03 


0.01 


0.01 


0.00 


0.00 


0.00 0.00 



TABLE 21: Coefficients of tfie nonperturbative parameters of tfie lepton energy OPE testing 
Moments, wfiere A = 5 [0.2, 1.3, 1, 1], D2 = ,5 [0.8, 1, 0.1, 1.3] , D3 = 5 [0.7, 1.6, 1.5, 1.5], 
D^^ = S [2.4, 1, 1.9, 0.8] and Dl^S [2.9, 1.4, 2.2, 1.3]. 





IS CKs^ a/ po 


as as^ l3o 


Combined Contributions 


pmin 


e Ai5 e elr.M 


e Ais e e%r^M 


e Ai5 e efj^^^ 


Di 
D2 

D3 


0.000 0.000 0.000 

-0.001 0.000 0.000 
0.000 0.000 0.000 
0.001 0.001 0.001 


0.000 0.000 0.004 
-0.001 -0.001 -0.012 
0.000 0.000 0.003 
-0.001 -0.001 0.003 


0.001 0.000 0.004 
-0.002 -0.001 -0.012 
0.000 0.000 0.003 
0.000 0.000 0.004 



TABLE 22: Coefficients of tfie perturbative parameters of tfie iepton energy OPE Moments. 



Labei 


Predicted Vaiue 


Measured Vaiue 


Di 


0.5459 ±0.0001 e ± 0.0010(N.P.) 




D2 


1.7585 ±0.006 e ± 0.0036(N.P.) 




D3 


0.5200 ±0.0001 e ± 0.0014(N.P.) 


0.5193 ± 0.0008 


Di 


0.6053 ±0.0002 e ±0.0018(N.P.) 


0.6036 ± 0.0006 



TABLE 23: Predictions and measurements for iepton energy OPE Testing Moments. 



[1] M. A. Shifman and M. B. Voloshin, Sov. J. Nucl. Pfiys. 41, 120 (1985); J. Chay, H. Georgi 

and B. Grinstcin, Phys. Lett. B247, 399 (1990). 
[2] I. I. Bigi et al.,Fhys. Lett. B 293, 430 (1992) [(E) ibid. B297 477 (1993)]; I. I. Bigi, et al, 

Phys. Rev. Lett. 71, 496 (1993). 
[3] A. V. Manohar and M. B. Wise, Phys. Rev. D 49, 1310 (1994); B. Blok et al. Phys. Rev. D 

49, 3356 (1994) [(E) ibid. D50 3572 (1994)]; T. Mannel, Nucl. Phys. B 413 396 (1994). 
4] M.B. Voloshin, Phys. Rev. D 51, 4934 (1995). 
5] M. Gremm and I. Stewart, Phys. Rev. D 55, 1226 (1997). 

6] M. Gremm, A. Kapustin, Z. Ligeti and M. B. Wise, Phys. Rev. Lett. 77, 20 (1996). 

7] A.F. Falk, M. Luke, M.J. Savage, Phys. Rev. D 53, 2491 (1996). 

8] A.F. Falk, M. Luke, M.J. Savage, Phys. Rev. D 53, 6316 (1996). 

9] M. Gremm and A. Kapustin, Phys. Rev. D 55, 6924 (1997). 
10] Adam F Falk, Nucl. Phys. Proc. Suppl. Ill, 3-13, (2002). 
11] S. Chen et al. (CLEO Collaboration), Phys. Rev. Lett. 87, 251807 (2001). 
12] D. Cronin-Hennessy et al. (CLEO Collaboration), Phys. Rev. Lett 87, 251808 (2001). 



28 



R. Briere et al. (CLEO Collaboration), hep-ex/0209024' 

B. Aubert et al. (BABAR Collaboration), hep-ex/ 02070841 
DELPHI Collaboration, ICHEP 2002, 2002-071-CONF-605. 
DELPHI Collaboration, ICHEP 2002, 2002-070-CONF-604. 

C. W. Bauer, Z. Ligeti, M. Luke, A.V. Manohar, Phys. Rev. D67, 054012 (2003). 
M. Battiglia et al., Phys. Lett. B556, 41-49 (2003). 

I. Gambino, N. Uraltsev, hep-ph 0401063 , 

C. W. Bauer, M. Trott, Phys. Rev. D 67,014021 (2003). 

Aida X. El-Khadra and M. Luke, Ann. Rev. Nucl. Part. Sci. 52:201-251 (2002). 

A. Hoang, Z. Ligeti and A.V. Manohar, Phys. Rev. D 59, 074017, (1999). 

A. Hoang, Z. Ligeti and A.V. Manohar, Phys. Rev. Lett. 82, 277, (1999); A. H. Hoang and 

T. Teubner, Phys. Rev. D 60, 114027 (1999) . 

I. I. Bigi, M. A. Shifman, N. G. Uraltsev and A. I. Vainshtein, Phys. Rev. D 50, 2234 (1994). 

M. Beneke and V. M. Braun,Nucl. Phys. B 426, 301 (1994). 

M. E. Luke, A. V. Manohar and M. J. Savage, Phys. Rev. D 51, 4924 (1995). 

M. Jezabek, J.H.. Kuhn, Nucl. Phys. B 320,20-44 (1989). 

G. Corbo, Nucl. Phys. B 212,99-108 (1983). 

A.F. Falk, M. Luke, M.J. Savage, Phys. Rev. D 53, 2491 (1996). 

F. De Fazio and M. Neubert, JHEP06, 017 (1999). 

A. Ah and E. Pietarinen, Nucl. Phys. B 154, 519, (1979). 

Y. Mr, Phys. Lett. B 221, 184 (1989). 

A. F. Falk, M. Luke, Phys. Rev. D 57, 424 (1998). 

M. Luke, M.J. Savage and M.B. Wise, Phys. Lett. B 343, 329 (1995). 

B. H Smith, M.B. Voloshin, Phys. Lett. B 340, 176 (1994). 
I. I. Bigi et al, Phys. Rev. D 52, 196, (1995). 
Valoshin, MB. Nucl. Phys. B 154, 365 (1979). 
Leutwyler. H, Phys. Lett. B 98, 447 (1981). 

Pineda. A, Nucl. Phys. B 494, 213 (1997). 

Vainstein Al, Zakharov VI, Shifman MA. JETP. Lett. 27, 55 (1978); Shifman MA., et al. 

Phys. Lett. B 77, 80 (1978); Shifman MA., et al. Phys. Lett. B 147, 385 (1979); 

Brambilla. N, Sumino Y, Vairo A. Phys. Lett. B 513, 381 (2001). 

M. Beneke and A. Singer, Phys. Lett. B 471, 233 (1999). 

A.Le Yaouanc et al., Phys. Lett. B488, 153 (2000), 

A.Le Yaouanc et al., Phys. Rev. D62, 074007 (2000), 

A.Le Yaouanc et al., Phys. Lett. B517, 135 (2001), 

I.I. Bigi and N. Uraltzev, Int. J. Mod. Phys. A16, 5201 (2001). 

D. Denson, I.I. Bigi and T. Mannel, Nucl. Phys. B 665, 367 (2003). 
N. Isgur, Phys. Lett. B 448, 111 (1999). 

A. Mahmo od et al. (CLEO Collaboration), Phys. Rev. D 67, 072001 (2003). 
M. Artuso, |hep-ph/0312270i 



29 



